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Abstract 

In this paper we implement Schwinger-Keldysh closed-time path integral formalism in non- 
equilibrium QCD to the definition of Collins-Soper fragmentation function. We consider a high pT 
parton in QCD medium at initial time tq with arbitrary non-equilibrium (non-isotropic) distribu- 
tion function f{p) fragmenting to hadron. We formulate parton to hadron fragmentation function 
in non-equilibrium QCD in the light-cone quantization formalism. It may be possible to include 
final state interactions with the medium via modification of the Wilson lines in this definition of 
the non-equilibrium fragmentation function. This may be relevant to study hadron production 
from quark- gluon plasma at RHIC and LHC. 
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I. INTRODUCTION 



RHIC and LHC heavy-ion colliders are the best facilities to study quark-gluon plasma 
in the laboratory. Since two nuclei travel almost at speed of light, the QCD matter formed 
at RHIC and LHC may be in non-equilibrium. In order to make meaningful comparison 
of the theory with the experimental data on hadron production, it may be necessary to 
study nonequilibrium-nonperturbative QCD at RHIC and LHC. This, however, is a difficult 
problem. 

Non-equilibrium quantum field theory can be studied by using Schwinger-Keldysh closed- 
time path (CTP) formalism l], [2I, ^. However, implementing CTP in non-equilibrium at 
RHIC and LHC is a very difficult problem, especially due to the presence of gluons in 



non-equilibrium and hadronization etc. Recent 



one-loop resumed gluon propagator in 



a. 



non-equilibrium in covariant gauge is derived in 

High pt hadron production at high energy e^e~, ep and pp colliders is studied by using 
CoUins-Soper fragmentation function . For a high pt parton fragmenting to hadron, 
CoUins-Soper derived an expression for the fragmentation function based on field theory and 
factorization properties in QCD at high energy . This fragmentation function is universal 
in the sense that, once its value is determined from one experiment it explains the data at 
other experiments. 

The derivation of parton to hadron fragmentation function in QCD medium based on first 
principle calculation is not done so far. This can be relevant at RHIC and LHC heavy-ion 
colliders to study hadron production from quark-gluon plasma. Further complication arises 
because the partons at RHIC and LHC may be in non-equilibrium. 

In this paper we note that, one can implement closed-time path integral formalism in 
non-equilibrium QCD to the definition of Collins-Soper fragmentation function. We con- 
sider a high Pt parton in medium at initial time Tq with arbitrary non-equilibrium (non- 
isotropic) distribution function f{p) fragmenting to hadron. We formulate parton to hadron 
fragmentation function in non-equilibrium QCD in the light-cone quantization formalism. 
The special case f(p) = -^^ — corresponds to the finite temperature QCD in equilibrium. 
This fragmentation function may be relevant to study hadron production from quark-gluon 
plasma at RHIC and LHC. 

We find the following definition of the parton to hadron fragmentation function in non- 
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equilibrium QCD by using closed-time path integral formalism. For a quark (g) with ar- 
bitrary non-equilibrium distribution function fq{k) at initial time, the quark to hadron 
fragmentation function is given by 



2z [1 + U{k)] J (2vr 



^troirac ^tr^oiorb^ < in\ij{x ,xt) ,xt] a^H{P+,OT)aH{P^,OT) $[0] ^(0)|m >] (1) 

where z i=^) is the longitudinal momentum fraction of the hadron with respect to the 
parton and Pt is the transverse momentum of the hadron. |m > is the initial state of the 
non-equilibrium QCD medium in the Schwinger-Keldysh in — in closed-time path formalism. 

For a gluon (g) with arbitrary non-equilibrium distribution function fg{k) at initial time, 
the gluon to hadron fragmentation function is given by 

1 f , - d'^~'^XT 

/ dx 

2zk+ [1 + fg{k)] J (27r) 

W[< ^n\F+'^ix-,XT) <f[x-,XT] aL(P+,OT)aH(P+,OT) $[0] F+(0)|2n >]. (2) 

° a=l 

The path ordered exponential 

POO 

^x%, = V exp[i^ / dXn- A^x^' + n^X) T^,] (3) 

is the Wilson line Qj. It may be possible to include final state interactions with the medium 

via modification of the Wilson lines in this definition of the non-equilibrium fragmentation 

□ 

function, similar to the pt distribution of the parton distribution function studied in [10l |. 

Eqs. ([T]) and O) can be compared with the following definition of CoUins-Soper fragmen- 
tation function |6|: 

D, 



' ' 2zJ (27r)'^-i 



'H/qyZ,r'T 

^troirac ^trcoiorb"" < 0|V'(x-,xt) $[x-,xt] a^, (P+ , OT)aH Or) $[0] 7/'(0)|0 >] (4) 



and 

^H/g{z,^T} 2zk+J {2ny-^ 

^E[<0|F+'^(x",xt) <^[x-,xt] aL(P+,OT)aH(P+,OT) $[0] P+(0)|0 >]. (5) 

° a=l 
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We will present derivations of eqs. ([T]) and ([2]) in this paper. 

The paper is organized as follows. In section II we briefly review the definition of CoUins- 
Soper fragmentation function in vacuum which is widely used at pp, ep and e^e~ colliders. 
In section III we give a brief description of Schwinger-Keldysh closed-time path integral 
formalism in non-equilibrium QCD. We implement closed-time path integral formalism in 
non-equilibrium QCD to Collins-Soper fragmentation function in section IV. Section V con- 
tains conclusion. 



II. COLLINS-SOPER FRAGMENTATION FUNCTION IN VACUUM 

In this section we will briefly review Collins-Soper fragmentation function in vacuum 
which is widely used at pp, ep and e~^e~ colliders. Consider a scalar gluon, for example, 
with four momentum k'^ in vacuum fragmenting to a hadron with four momentum P^. For 
application to collider experiments it is convenient to use light cone quantization formalism. 
The scalar gluon field 0(a;) can be written as 

1 r rlk~^ 

where and a are the creation and annihilation operators respectively, d = A — 2t where 
3 — 2e is the space dimension. The single particle parton state is given by 

\k^,kT>= a\k+,kT)\Q>, with a(A;+, A;t)|0 >= (7) 

with the normalization 

< krW^, k'r >= {2nY-^6{k+ - k'^)6^~\kT - k'r). (8) 

Note that the correct interpretation of the state |A; > is created by an appropriate Fourier 
transform of the corresponding field operator and should not be associated with on-shell 



condition A;^ = of the massless quark or gluon 12|. 

Consider the inclusive production of a hadron H created in the o-ut— state \H,X > from 
the parton a in the m— state \k > with the probability amplitude 

<H,X\k>. (9) 
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The probability distribution hk{P) of the hadron with momentum P from the parton of 
momentum k can be found from the above amphtude. Exphcitly 

hk{P) < k\k' >=J2 < k\H,X >< H,X\k' >=J2 < k\a^H{P)\X >< X\aH{P)\k' > 



X X 



=< k\a}H{P)(^H{P)\k' > (10) 

where a}^ is the creation operator of a hadron H. In the hght-cone quantization formahsm 
we find (by using eqs. (171) and ([H])) 

h{z, Pt) < k+, krlk'^, k'r >= 2z{2txY-^ DH/a{z, PT){2n)'^-^6{k+ - k'+)6'^-\kT - k'r) 
=< 0\a{k~^,kT)aUP'^,PT)aH{P^,PT)a\k'^,k'T)\0 > (11) 

where Dn/aiz, Pt) is the fragmentation function and z = is the longitudinal momentum 
fraction of hadron H with respect to parton a. Using 

1 r dk^ r dk'^ 

(27,)'-' < 0\<P{x-,XT)aUP\PT)aH{P\PT)m\0 >= J ^y^^'-'kr J ^y^d'-'k'. 



[< 0\e-'''-'^a{k+,kT)a^HiP^^PT)aH{P^,PT)aHk'^,k'T)\0 > 



x+=0 



(12) 



we find from eq. ( jTT|) 

DH/a{z,PT) = ^/ dx-^^^e^'^''^--"'^-^^ < 0\(j){x-,XT)aUP^,PT)aH{P+,PT)m\0 > ■ (13) 

It is convenient to rewrite the definition in a form analogous to the definition of the distri- 
bution of partons in a hadron. The transverse momentum is of the parton relative to the 
hadron rather than vice versa. For this purpose we make a Lorentz transformation to a 
frame where the hadron's transverse momentum is zero: 

{P\PT)^iP\0) 

(fc+,0)-^(fc+,-PTA). (14) 



Hence we find from eq. (1131) 

k^ f d'^~'^XT 



DH/a{z,PT) = —J dx-j^^;^e'^^^-+'''^-^^'' < 1 0(x- , XT)aL (P+, OT)aH (P+ , Ot) 0(0) 1 > . (15) 
The Pt integrated fragmentation function is given by 

dH/a{z) = J d^-^PTDH/giz.PT) 
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J rfx-e^^^^"/"[< 0\(j){x-)a^HiP^, OT)aH{P^, Or)0(O)|O >]. (16) 



A. Quark Fragmentation Function 



Following similar steps as above but performing calculation for quark we find the quark 
fragmentation function 



ttY-^ J 



itrD„„ itr„,„,|7+ < 0\<li(x-,XTW„(P+,OT)a„(P+,OTym\0 >] (17) 
where ip{x) is the quark field. The Pt integrated quark fragmentation function becomes 
dndz) = ^ / rfx-e^^+^^/^itroirac ^tr,oior[7+ < 0|V^(x-)al,(P+, Or)aH(P+, Ot)^(0)|0 >].(18) 



B. Gluon Fragmentation Function 

Following the above steps but for gluons we find the gluon fragmentation function 

^E[< 0|F+'^(x-,a:T)akP+,OT)aH(P+,OT)F+(0)|0 >] (19) 

° a=l 

where 

F^" = df'Al - d^Ai^. (20) 
The Pt integrated gluon fragmentation function becomes 

dH/,{z) = / dx-e^^^^-/^j:[< 0\F^^^{x~)aUP\OT)an{P\OT)F^M\0 >]■ (21) 

a=l 

C. Wilson Lines and Fragmentation Functions 

The quark and gluon fragmentation functions as defined above are not gauge invariant. 
Gauge invariant parton fragmentation functions are obtained by incorporating Wilson lines 
$[x]afc into the definition of the quark and gluon fragmentation functions [sl, 1^, 1^, 11 1. We 
find 

DH/g{z,PT) = ^^^^^ I dx-d'^-'xre^'^^-^'''^-^^/' 

^troirac ^treoiorb^ < 0|V^(x-,xt) ^x- , xt] al,(P+, OT)a//(P+, Ot) $[0] ^(0)|0 >] (22) 
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where is given by eq. ([3]) with T"'' in the fundamental representation of SU(3). 

Similarly, incorporating Wilson lines, we find the gauge invariant gluon fragmentation 
function 

i^[< 0\F+^{x-,xt) aUP^,OT)aH{P^,OT) $[0] F+(0)|0 >] (23) 

° a=l 

where is given by eq. ([3]) with T"'' = /"^^^ in the adjoint representation of SU(3). 



III. NON-EQUILIBRIUM QCD USING CLOSED-TIME PATH FORMALISM 

Unlike pp collisions, the ground state at RHIC and LHC heavy-ion collisions (due to the 
presence of a QCD medium at initial time t = tin (say tm=0) is not a vacuum state |0 > 
any more. We denote \in > as the initial state of the non-equilibrium QCD medium at tin- 
The non-equilibrium distribution function f{k) of a parton (quark or gluon), corresponding 
to such initial state is given by 

< a\k)a{k') >=< tn\a\k)a{k')\tn >= f{k){2TiY-^5^'^-^\k-k') (24) 

where we have assumed space translational invariance at initial time. 

Finite temperature field theory formulation is a special case of this when f{k) 
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Consider scalar gluons first. In the CTP formalism in non-equilibrium there are four 
Green's functions 

G++(x,x') =< in\T(j){x)(j){x')\in >=< T0(x)0(x') > 

G {x,x') =< in\T(j){x)(j){x')\in >=< T(f){x)(j){x') > 

G-\ — {x,x') =< in\(f){x')(f){x)\in >=< (f){x')(f){x) > 

G — \-{x,x') =< in\(j){x)(j){x')\in >=< 0(x)0(x') > (25) 



where +(— Isign corresponds to upper(lower) time branch of the Schwinger-Keldysh closed- 
time path y]- T is the time order product and T is the anti-time order product. The 
field 0(x) is in Heisenberg representation. Explicitly 

T(j){x)(j){x') = e{t - t')(j){x)(i){x') + e{t' - t)(j){x')(j){x) 

f(j){x)(f){x') = e{t' - t)(j){x)(f){x') + e{t - t')(f){x')(f){x). (26) 



At initial time t = tin = the Heisenberg picture coincide with the Schrodinger and 
interaction pictures. We write 

where a\k) and a{k) are creation and annihilation operators respectively. The commutation 
relations are given by 

[a{k),a\k%=o = {27rY-'S^''-'\k-k'), 

[a{k),a{k%=o = [at(fc),at(fc')]t=o = 0. (28) 
We assume space-translational invariance at initial time and find 

[G,j{x, x%=o = [ / d''kG,,{k)e-'''<''--\=o (29) 
where i,j are +, — . Explicitly 

d'^^^k r d'^~^k' 

G_+{x,x') =< (j){x)(j){x') >-- 



/2F(27r)^-i^ V2F^(27r)'^-i 
< tn\[a{k)e-'^-'' + a^(fc)e'^-"] [a(fc')e"''''-"' + a\k')e"''-''']\m > . (30) 

Using Bogolyubov transformation we can set 

< in\a{k)a{k')\in >=< in\a'^ {k)a'' {k')\in >= 0. (31) 

By using eqs. (I21D, (EH]) and ([SI]) in eq. ^ we find 

[G^^ix,x')U = [J ^J12y-^ + /(fc)]e-^-(--') + /(fc)e-^-(--')],=o 

"^'^ -5(fc2)e-^^ (^-^') [9{ko) [1 + fm + e{-ko)f{-k)U = [ / d'kG^+{k)e-^'<^-^\=o. 



{2tt 



(32) 



Similarly 



[G+4k)]t=o = S{e)[9{~ko) + 9{ko)f{k) + 9{-ko)f{-k)] 
[^++(^)]*=o = + ^ik')mo)f{k) + e{-ko)f{-k)] 

[G-ik)]t=o = + 6iemk,)fik) + e{-k,)f{~k)]. (33) 



A. Quarks in non-Equilibrium 



The non-equilibrium (massless) quark propagator at initial time t = tin is given by 
(suppression of color indices are understod) 



-27i5ik^)e{-ko) + 27iS{k^)fgik) 
h 2n6{e)f,{k) 



(34) 



fc2- 



where where i,j = +, — and fq{k) is the arbitrary non-equilibrium distribution function of 
quark. 



B. Gluons in Non-Equilibrium 



We work in the frozen ghost formalism 



B,|5 



5| where the non-equilibrium gluon propagator 

at initial time t = tin is given by (the suppression of color indices are understood) 

k^^k^ 

C^^ik),, = -zig'^'^ + {a- 1)-^] GTf (A;) - zT^'^G^'^^ik) (35) 

where i,j = +, — . The transverse tensor is given by 

{k-u){u^k'^ + u'^k^)-k^k--k'u^u'^ 
^ ^^'-^ {k ■ uf - P ^"^^^ 

with the flow velocity of the medium u'^. Gi^{k) are the usual vacuum propagators of the 
gluon 



Girw 



_1_ -2n5ie)ei-ko) 

2\Q{U \ 1 



\-2nSik')9iko; -p- 
and the medium part of the propagators are given by 



(37) 



G^'^^k) = 2n6ie)f,ik) I I I I . (38) 



IV. FRAGMENTATION FUNCTION IN NON-EQUILIBRIUM QCD USING 
CLOSED-TIME PATH FORMALISM 

For simplicity, let us consider the scalar gluon first in the light cone quantization formal- 
ism. Generalizing the vacuum analysis in we define the state \k~^,kT > of the fragmenting 
gluon in non-equilibrium QCD medium at initial time x"*" = xf^ (say at x^^ = 0) 

\k+ , kT >= {k^ , kT)\in > . (39) 
9 



The non-equilibrium distribution function f{k~^, ^t) of the fragmenting gluon is given by 

< in\a\k+, kT)a{k'^, A;^)|m >= f{k+, A;r)(2vr)^-i(5(A;+ - k'^)6^'^-^\kT - k'^) (40) 

where we have assumed the space {x~ and xt) translational invariance at initial time x'^ = 
xfn- The commutation relations are given by 

[a{k+,kT),a\k'^,k'j,)l+=o = {2ny-'6{k+ - k'^)6^''-^\kT - k'^), 

[a{k+, kr), a{k'+, A;^)].+=o = Wik+, kr), a\k'^ , A;^)].+=o = (41) 

By using eqs. (!24l) and (!28|) we find 

< k~^, kT\k'~^ , k'x >=< in\a{k~^, kT)a\k'^ , k'T)\in > 

= i2nf-'6{k+ - k'^)6'-\kT - k'T)[l + f{k+,kT)]. (42) 

Consider the inclusive production of hadron H created in the owt— state \H,X > from a 
scalar gluon in non-equilibrium in the initial state \k > with the probability amplitude 

<H,X\k>. (43) 

X being other outgoing final state particles. Similar to the vacuum case of Collins-Soper 
fragmentation function, the correct interpretation of the above state |fc > is created by an 
appropriate Fourier transform of the corresponding field operator and should not be asso- 
ciated with on-shell condition A;^ = of the massless quark or gluon [l^. The distribution 
hk{P) of the hadron H with momentum P from the parton of momentum k can be found 
from the above amplitude. We find 

^ < A;, kT\H, X >< H,X\k-^', k'^ >= hk{P) < k+, kT\k+\ k'j, > ■ (44) 

X 

For the left hand side we write 

< k+,kT\H,X >< H,X\k+\kT>=Y. <k+,kT\a^H{P)\X >< X\aH{P)\k+',k'T> 

X X 

=< k+, kT\a^HiP)(^H{P)\k^\ k'j, > . (45) 

Equating eqs. and fH51) and by using eq. (15^ we find 

< in\a{k+,kT)a^HiP^^PT)aH{P^,PT)a^{k'^,k'T)\in > 

= 2z{2ny-^DH/a{z, Pt) < in\a{k+, kT)a\k'^ , k'r) \in > . (46) 
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This expression is exactly similar to that of the jet fragmentation function in vacuum as 
given in eq. ffTTl) except that the vacuum expectation is replaced by medium average at the 
initial time x"*" = x^. Using eqs. (15^ and (HUi) we find 

< in\a{k+, kT)a^H{P^, PT)aH{P^, PT)a\k'^, >= 2z{2tiY-^ Du/aiz, Pt) [1 + /(A;+, M]- 

(47) 

From eq. (1271) we obtain 

i2nY-' < ^n\<P{x~,XT)aUP\PT)aHiP\PT)m\^n >= J -^d'-'^T J -^d'-'k'r 

[< in\e-'''-''a{k+,kT)a^HiP^,PT)aH{P^,PT)a\k'^,k'T)\in >]^+=o. 



< tn\(j)a{x-,XT)aUP^,PT)aH{P^,PT)<l)am^n > . 



From eq. ([T4j) we find 



Using this in eq. f H7j) we find the expression of the fragmentation function in non- 
equilibrium QCD 

f d'^~'^XT Jk+X--ikr, 



(49) 



'^/'^^"'"^^ = z[l + f{k^,kT)]J 

< m|0a(x-,xr)al^(P+,OT)aH(P"*',OT)0a(O)|m > . 

(50) 

In the above expression, f{k~^, kx) is the non-equilibrium distribution function of the frag- 
menting gluon at initial time = xf^ and \in > is the initial state of the non-equilibrium 
QCD medium in the Schwinger-Keldysh in — in closed-time path formalism. 

A. Quark Fragmentation Function in Non-Equilibrium 

Following the above steps, but for quarks, we find the quark fragmentation function in 
non-equilibrium QCD 

( p ) ^ I f ^ - d'^'^XT ik+x-+iPr-xr/z 

' 2z [1 + /,(fc+, M] i i^^r-' 

2trDirac gtrcoiorb"^ < m| ?/'(x" , XT)aJ;^ (P+, Ot)^^ (^'^ , Ot)^ (0) |m > (51) 



(48) 
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where fq{k~^, kx) is the non-equihbrium distribution function of the fragmenting quark at 
initial time. 



B. Gluon Fragmentation Function in Non-Equilibrium 

For gluons we consider frozen ghost formahsm described above (135|1 . Hence all the analysis 
of the above can be applied. Carrying out the similar algebra as above we find the gluon 
fragmentation function 

Dul (Z Pt) = ^ / dx- ik+x-+iPT-XT/z 

W[< m|F+'^(x-,XT)aUP+,OT)a^(P+,OT)F+(0)|m >] (52) 



a=l 



where fg{k~^,kT) is the non-equilibrium distribution function of the fragmenting gluon at 
initial time. 

C. Wilson Lines 

The above parton to hadron fragmentation function definition in non-equilibrium QCD 
is not gauge invariant. To make it gauge invariant we need to incorporate Wilson lines. 
Incorporating the Wilson lines (eq. ([3])) into the definition of the fragmentation function eq. 
( ISTl) we find the gauge invariant quark fragmentation function in non-equilibrium QCD 



^troirac ^treoiorb^ < m|V^(x-,xr) <^[x-,xt] al^(P+,OT)aH(P+,Or) $[0] ^(0)|m >] 

(53) 

where fq{k~^, kx) is the non-equilibrium distribution function of the fragmenting quark at 
initial time x^ = xf^. \in > is the initial state of the non-equilibrium QCD medium in the 
Schwinger-Keldysh in — in closed-time path formalism. This reproduces eq. ([T]). 

Incorporating the Wilson lines (eq. ([3])) into the definition of the fragmentation function 
eq. ( l52l) we find the gauge invariant gluon fragmentation function in non-equilibrium QCD 
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\Y.\<MFa'i^~.^T) $[x-,xt] al,(P+,OT)aj,(P+,OT) <f[0] F+(0)Kn >] (54) 



a=l 



where fg{k~^,kT) is the non-equihbrium distribution function of the fragmenting gluon at 
initial time = xf^. \in > is the initial state of the non-equilibrium QCD medium in the 
Schwinger-Keldysh in — in closed-time path formalism. This reproduces eq. ([2]). 

It may be possible to include final state interactions with the medium via modification 
of the Wilson lines in this definition of the non-equilibrium fragmentation function, similar 



to the pt distribution of the parton distribution function studied in 



V. CONCLUSIONS 



3- 



In this paper we have implemented closed-time path integral formalism in non-equilibrium 
QCD to the definition of Collins-Soper fragmentation function. We have considered a high 
Pt parton in QCD medium at initial time tq with arbitrary non-equilibrium (non-isotropic) 
distribution function f{p) fragmenting to hadron. We have formulated parton to hadron 
fragmentation function in non-equilibrium QCD in the light-cone quantization formalism. 
This may be relevant to study hadron production from quark-gluon plasma at RHIC and 
LHC. It may be possible to include final state interactions with the medium via modification 
of the Wilson lines in this definition of the non-equilibrium fragmentation function, similar 
to the Pt distribution of the parton distribution function studied in [10|]. This will be the 
subject of a future analysis. 
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